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Abstract 

We study static spherically symmetric dyonic black holes in Einstein- Yang- 
Mills-Higgs theory. As for the magnetic non-abelian black holes, the domain of 
existence of the dyonic non-abelian black holes is limited with respect to the 
horizon radius and the dimensionless coupling constant a, which is proportional 
to the ratio of vector meson mass and Planck mass. At a certain critical value 
of this coupling constant, a, the maximal horizon radius is attained. We derive 
analytically a relation between a and the charge of the black hole solutions and 
confirm this relation numerically. Besides the fundamental dyonic non-abelian 
black holes, we study radially excited dyonic non-abelian black holes and globally 
regular gravitating dyons. 
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1 Introduction 



SU(2) Yang-Mills-Higgs (YMH) theory, with the Higgs field in the adjoint represen- 
tation, possesses globally regular particle-like solutions, such as the 't Hooft-Polyakov 
magnetic monopole |]l| and the Julia-Zee dyon, which carries both magnetic and electric 
charge 0. 

In SU(2) Einstein- Yang-Mills-Higgs (EYMH) theory, for small gravitational con- 
stant the gravitating fundamental monopole 0, H, ^ and the gravitating fundamental 
dyon [0 emerge smoothly from the corresponding fiat space solutions. The mass of the 
gravitating fundamental monopole [§, |^ and dyon [^] solutions decreases with increas- 
ing gravitational constant, corresponding to increasing coupling constant a, which is 
proportional to the ratio of vector meson mass and Planck mass. The regular solutions 
cease to exist beyond some maximal value of a, depending on the electric charge Q 

Besides these globally regular solutions, in SU(2) EYMH theory magnetic [1, i|, H, 
and dyonic non-abelian black hole solutions exist. They emerge from the globally 
regular solutions, when a finite regular event horizon Xh is imposed. Consequently, they 
have been characterized as "black holes within magnetic monopoles" |^ or "black holes 
within dyons" . Distinct from the corresponding embedded Reissner-Nordstr0m (RN) 
black holes |TD|, these non-abelian black hole solutions represent counterexamples to 



the "no-hair" conjecture, because they carry non-trivial non-Coulomb-like fields outside 
their horizon. In contrast, pure SU(2) EYM theory possesses neither regular dyon 
solutions nor dyonic black holes other than embedded RN solutions [|12|, |13|. 



The magnetic non-abelian black holes exist only in a limited domain of the a — 
Xh plane |Q, ^ ^. The domain of existence consists of two regimes, where for fixed 
coupling constant a and varying horizon radius x/j, the solutions approach limiting 
solutions in two distinct ways. At the transition point between these two regimes, 
given by the critical value d, the maximal horizon radius is attained 0, ^. Limiting 
our investigations to vanishing Higgs self-coupling, we here determine the domain of 
existence of the dyonic non-abelian black holes, finding also two regimes for the coupling 
constant a. 

Besides the fundamental gravitating monopole and dyon solutions there are radi- 
ally excited solutions ||, H, where the gauge field function possesses nodes. Both 
the radially excited monopole and dyon solutions are related to the globally regular 
Einstein- Yang-Mills (EYM) solutions, found by Bartnik and McKinnon [14|, and, like 
these solutions, have no fiat space counterparts. Besides the radially excited globally 
regular dyon solutions, we here consider the radially excited dyonic non-abelian black 
hole solutions and determine their domain of existence. 

In section II we derive the EYMH equations of motion. We present the equations 
both in Schwarzschild-like coordinates, employed in the numerical calculations , and 
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in isotropic coordinates, to extend the analytical calculations 0] to the dyonic case. In 
section III we briefly review the embedded RN solutions. We discuss the globally regular 
dyon solutions in section IV. In addition to vanishing Higgs self-coupling, we consider 
the fundamental solutions also for finite Higgs self-coupling, in order to demonstrate 



that, as in flat space |T^, a maximally possible charge Qmax arises. In section V we 
discuss the dyonic non-abelian black hole solutions. In particular, we derive the critical 
value a for the dyonic solutions analytically and confirm it numerically. For vanishing 
Higgs self-coupling, we determine the domain of existence for the fundamental and the 
first radially excited dyonic non-abelian black hole solutions, and we discuss the pattern 
how the black hole solutions approach limiting solutions. We conclude in section VI. 
In Appendix A, finally we show the local existence of the non-abelian regular dyon and 
dyonic black hole solutions at the origin and at the horizon, respectively. 



2 Einstein- Yang-Mills-Higgs Equations of Motion 

We consider the SU(2) EYMH action 

S = Sg^ Sm = J LcV^d^x + J LuV^d^x (1) 

1 



with 



(2) 



and 



where 



Lm = -\f;,f'^^'' - \D,rD^r - ^{rr - , (3) 

F^u = d,K - QuAl + ge-'^A\Al , (4) 

/^^r = + S^^'^Alc^' , (5) 

g is the gauge coupling constant, /3 is the dimensionless Higgs coupling constant, pro- 
portional to the ratio of Higgs boson mass and vector boson mass, and v is the Higgs 
field vacuum expectation value. Variation of the action eq. (|^) with respect to the 
metric (yf^jy, the gauge field A"-^ and the Higgs field 0" leads to the Einstein equations 
and the matter field equations. 

2.1 Ansatz 

For static spherically symmetric globally regular and black hole solutions the metric 



can be parametrized as H, 16 



ds^ = g^.dx^dx" = -e^^^'^^dt^ + e^^^'^UR^ + r\R){de^ + sin^ ed(P^) . (6) 
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For the gauge and Higgs field we employ the spherically symmetric ansatz 

At = erJ{r)v , (7) 

1 - Kir) - 1 - Kir) 
A = , Afi = -e^ ^ , A^, = eft — sin 6 , (8) 
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and 



erHi^r)v , 



(9) 



with the standard unit vectors e^, and e^. For A" = 0, the solutions carry only 
magnetic charge 0, P]. 

With these ansatze we obtain the action 



S = ^ I dtdRe^^^ 



1 + e-'^ [{rY + v'{r') 



+47rG 



2„,2 



-2A 



\K'f r2 



2 +-mvi 



(10) 



where z/, A, r, i^', J and H are functions of i?, the prime indicates the derivative with 
respect to R, and 



[l-K 

2^2j,2 



2\2 



I t^2tt2 2 , 9 2/tt2 i n2 4 

+ K H V H r [H — 1) v 



'IV. 



2.2 Schwarzschild-like coordinates 

We first employ Schwarzschild-like coordinates, corresponding to the gauge r{R) = R. 
Renaming R = r, the spherically symmetric metric reads 



ds'^ = -A^Ndt^ + N-^dr"^ + r'^{de'^ + sin^ Odcf)'^) 
with the metric functions 



A{r) = e 
Nir) = 1 _ MO ^ ^-2A 



(12) 

(13) 
(14) 



and the mass function m{r). 

We now introduce the dimensionless coordinate x and the dimensionless mass func- 
tion n, 

X = gvr , /i = gvm , (15) 
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as well as the dimensionless coupling constant a 

0? = AttGv^ . (16) 

The tt and rr components of the Einstein equations then yield the equations for 
the metric functions, 

+ + H'K' + ^x\H' - 1)A , (17) 



2x^ 4 



and 



A' = o?x , H — + M 



I ^2JV22,2 3,2 



where the prime now indicates the derivative with respect to x. For the matter functions 
we obtain the equations 

(AiVir')' = Air(^ + i/^-^) , (19) 



'x'J'X 2JK^ 



(20) 



\ A J AN ' 
and 

{x^ANHj = AH (2K^ + P^x^H'' -1)) . (21) 
Introducing the new function 

P=j (22) 

and using eq. (|18|), we see that the metric function A can be eliminated from the set of 
equations ([T7|) , (|T9|)-(pl|). These can then be solved separately, with the metric function 
A being determined subsequently by (|TB|). 

In the following we refer to the above set of functions, variables and parameters, 
also employed in H], as notation I. 



2.3 Isotropic coordinates 

To make contact with the analytical results of P, || and to extend them to the case 
of gravitating dyons and dyonic black holes, we now present the equations of motion 
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in the notation of [Q, in the following referred to as notation II. In the Schwarzschild 
gauge the following identifications hold 

[6] : a; A N K H J a 

§ : ar e^+^ e'^^ W aH J a ^ ^ ' 

Along with 0, we now employ the gauge choice = corresponding to isotropic 
coordinates, and introduce the coordinate 

r = lni?, (24) 

and the dot ( ' ) now indicates the derivative with respect to r. To obtain a system of 
first order equations we further employ for the first derivatives the new variables ^ 

r W 

N= - , K = iy + N , U= — , V = H , (25) 

r r 

and we introduce 

B = e-^J , C = B . (26) 

(Note, that the symbol here has a totally different meaning from eq. (0).) 
Then the following system of first order autonomous equations is obtained 



r 


= Nr 


(27) 


N 


= N{k -N)- (2f/2 + 1/2) _ 2a'^BV 


(28) 


k 


= 1 K^ + 2U^ ^^''\h^ ay 2W'H^ + 2a^W'B' 


(29) 


H 


= V 


(30) 


V 


= ^^''\h^ a^)H + 2W^H kV 


(31) 


W 


= rU 


(32) 


U 


r 


(33) 


B 


= C 


(34) 


C 


= -kC + B{2W'^ + 






+2W^H^ + ^'"'^ {H^ a'f) Aa'B^W'' 


(35) 



supplemented by the constraint 

2kA^ = l + N^ + 2U'^ + V'^ - 2V2 + 2a^W'^B'^ - a^{C + B{k - A^))^ (36) 
compatible with the equations. 



6 



3 Embedded Reissner-Nordstr0m Solutions 



As noted long ago [jTO|, the Einstein- Yang- Mills-Higgs system admits embedded RN 
solutions. In notation I, RN solutions with mass /ioo, unit magnetic charge and arbitrary 
electric charge Q are given by 

= yUoo - ^ ^ ^ A{x) = 1 , (37) 

K{x) = , J(x) = J^~^, H{x) = 1 , (38) 

X 

where J^o is a priori arbitrary, but here chosen as Jqo = Q/xh (see eq. (0)). At the 
regular horizon Xh the metric function 

Nix) = -'-^-^^o^+/(^ + Q') (39) 

vanishes, N{xh) = 0, yielding 

Xh = fioo + v//i^-«2(i + g2) (40) 

or equivalently 

xl + a\l + Q^) 
= 2^, • 

For fixed values of a and Q, RN solutions exist for Xh > <y\/l + Q^, independent of /3. 
In the extremal case 



a;/, = /ioo = ttyi + , (42) 

extremal RN solutions are obtained. Extremal RN solutions are characterized by Q and 
a. In particular, for extremal RN solutions the function B{x) defined in (|2^) becomes 
constant ^ ^ 

^ ^ a7n ^ aVT+W ■ ^^^^ 



4 Dyon Solutions 

Let us first consider the globally regular particle-like solutions of the SU(2) EYMH 
system. 

In the Prasad-Sommerfield limit, (3 = 0, the dyon solutions in fiat space are known 
analytically P|, whereas for finite (3 they are obtained only numerically ll5|. In the 
presence of gravity, the corresponding gravitating monopole and dyon solutions are 
obtained only numerically as well. The gravitating dyon solutions have many features 
in common with the gravitating monopole solutions [0. 
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After presenting the boundary conditions for the asymptotically flat globally regular 
solutions, we here briefly discuss the fundamental dyon solutions both in the Prasad- 
Sommerfield limit and for finite /3. Then we turn to the excited dyon solutions. 

4.1 Boundary conditions 

Requiring asymptotically flat solutions implies in notation I that the metric functions 
A and n both approach a constant at infinity. We here adopt 

A{oo) = 1 , (44) 

and yu(oo) = fioo represents the dimensionless mass of the solutions. The matter func- 
tions also approach constants asymptotically, 

K{oo) = , J{oo) = Joo , ^(oo) = 1 , (45) 

where for magnetic monopole solutions Joo = 0. The asymptotic fall-off of the function 
J{x) determines the dimensionless electric charge Q (see eq. (^8])). 
Regularity of the solutions at the origin requires 

m = , (46) 

and § 

K{0) = 1 , J(0) = , H{0) = . (47) 

The local existence of a family of analytic solutions obeying these conditions is shown 
in Appendix A.l. 

4.2 Fundamental dyons 

/3 = 

Like the gravitating monopole solutions, the gravitating dyon solutions exist up to 
a maximal value of the coupling constant a. Beyond this value no dyon solutions exist. 
The fundamental dyon branch does not end at the maximal value Omax, but bends 
backwards and extends up to the critical value ac- Since = AirGv'^, variation of 
the coupling constant along the fundamental dyon branch can be considered as first 
increasing G with v fixed up to the maximal value amax? and then decreasing v with G 
fixed up to the critical value ac- 

Completely analogously to the fundamental monopole branch Q, the fundamental 
dyon branch reaches a limiting solution at the critical value ac, where it bifurcates with 
the branch of extremal RN solutions of unit magnetic charge and electric charge Q . 
The critical value a^Q) depends only slightly on the charge Q 
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For a — i> the minimum of the metric function N{x) of the dyon solutions decreases 
monotonically, approaching zero at Xc = dcV^ + Q^- For x > Xc, the metric function 
N{x) of the hmiting solution corresponds to the metric function of the extremal RN 
black hole with horizon 

VI + <5^, unit magnetic charge and electric charge 
Q. Likewise, for x > Xc, the other functions of the limiting solution correspond to those 
of this extremal RN black hole P] . Consequently, also the mass of the limiting solution 
coincides with the mass of this extremal RN black hole. 

On the interval < x < Xc the functions N{x), K{x) and H{x) of the limiting 
solution vary smoothly from their respective boundary values at x = to those of 
the extremal RN solution at Xc, whereas the function J{x) is identically zero. The 
metric function A{x) is identically zero, as well, on the interval < x < Xc, but 
discontinuous at Xc- Consequently, in the limiting spacetime the inner and outer parts 
are not analytically connected ^. 

/3 > 

Let us briefly consider the dyon solutions for finite Higgs self-coupling, (3 > 0, 
because we here encounter a phenomenon, also present for the black hole solutions at 
vanishing Higgs self-coupling, (3 = 0. 

In fiat space dyon solutions of arbitrarily large charge exist only for j3 = 0. In 
contrast, for finite values of (3 dyon solutions exist only up to a maximal value of the 
charge, Qmax(/3) [|15[. As the maximal value of the charge is approached, the function 



J(x) approaches asymptotically the limiting value 

lim = 1 • (45 



ii max 



At this point the solutions change character and become oscillating instead of asymp- 
totically decaying [p!5[] . 

As illustrated in Figs. 1-2, this phenomenon persists in curved space. In Fig. 1 we 
show the mass of the flat space dyon (solid line) and of the gravitating dyon (dashed 
line, corresponding to a = 0.5) as a function of the charge Q for = and = 0.5. 
In Fig. 2 the corresponding asymptotic values of the function J(x) are presented. We 
note, that for finite (3 also gravitating dyon solutions exist only up to a maximal value 
of the charge, where Jqo — >■ 1- 



4.3 Radially excited dyons 

Besides the branch of fundamental dyon solutions there are branches of radially excited 
dyon solutions, where the gauge field function K{x) of the n-th radially excited dyon 
solution has n nodes. 

These solutions have no fiat space counterparts, and the variation of a along a 
branch of excited solutions must be interpreted as a variation of v with G fixed. In 
the limit a —>■ the Higgs field vacuum expectation value therefore vanishes, while G 
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remains finite. Because of tlie particular clioice of dimensionless variables (]T5|), in this 
limit the solutions shrink to zero size and their mass /i diverges. As for the radially 
excited monopole solutions the coordinate transformation x = x/a leads to finite 
limiting solutions in the limit a ^ 0, which are the Bartnik-McKinnon solutions [Q. 

In Fig. 3 we show the normalized mass /ioo/« as a function of a for the first excited 
dyon branch with electric charge Q = 1 and /3 = 0. For comparison also the first 
excited monopole branch is shown. 

Like the radially excited monopole solutions, the radially excited dyon solutions 
exist only below a critical value of the coupling constant a. For the radially excited 
monopole solutions this critical value is Uc = v^/2 |^. For the radially excited dyon 
solutions the critical value is larger and will be discussed in the context of the dyonic 
black hole solutions. 



5 Black Hole Solutions 

We now turn to the dyonic black hole solutions of the SU(2) EYMH system, choosing 
/3 = 0. Again, the SU(2) EYMH dyonic black holes have many features in common 
with the magnetic non-abelian black holes In particular, dyonic non-abelian black 
hole solutions exist in a limited domain of the a-Xh plane, which depends on the charge 
Q, and, in the limit Xh — > 0, the corresponding globally regular solutions are obtained. 
But a new phenomenon occurs for small values of a, where the non-abelian solutions 
no longer bifurcate with non-extremal RN solutions. 

In the following, we briefly present the boundary conditions. We generalize the 
considerations of to dyonic solutions and derive analytically the critical value of the 
coupling constant, a. We then present our numerical calculations for the fundamental 
dyonic black hole solutions and their radial excitations. 



5.1 Boundary conditions 



Imposing the condition of asymptotic flatness, the black hole solutions satisfy the same 
boundary conditions at infinity as the regular solutions. The existence of a regular 
event horizon at Xh leads in notation I to the conditions for the metric functions 



and A{xh) < oo, and for the matter functions 



2 



N'K'\ 



K 



x" 



(49) 

(50) 
(51) 
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and 



x'N'H'l^ = H (2K^ + f3^x\H^ - 1) 



(52) 



The local existence of a family of solutions obeying these conditions at the horizon is 
demonstrated in Appendix A. 2. 



5.2 Critical value of a 

For magnetic black holes, there are two regimes of the coupling constant a, corre- 
sponding to two distinct patterns of reaching a limiting solution as a function of the 
horizon radius Xh- For larger values of a, the solutions bifurcate with an extremal RN 
black hole solution, when the horizon radius appoaches its critical value xj^. Here the 
solutions tend towards the corresponding extremal RN solution with horizon radius 
x^^ = a only on the interval x > x]^^, whereas they tend to a generic non-abelian 
limiting solution on the interval x'^ < x < x^^ In contrast for smaller values of 
a the solutions bifurcate with a non-extremal RN solution on the full interval x > 

1^ . The transition point between these two regimes corresponds to the critical value 
a, where the non-abelian black hole solutions tend towards the extremal RN solution 
with horizon radius x^^ = a on the full interval x > x'f^ = xj^^, i.e. x'f^ and xf^^ coin- 
cide. For the magnetic non-abelian black holes this critical value is given by d = 

Analogously, for dyonic non-abelian black holes there are also two regimes of the 
coupling constant a, with a marking the transition point between these two regimes. 
We now derive this critical value a for the dyonic non-abelian black hole solutions, 
following closely the reasoning of M. 



Let us then consider the fixed points of the system (|27|)-(|35D. Employing notation 
II, the relevant fixed point is given by the configuration 

Hh = a , Wh = , K/i = 1 , rh= , , (53) 

^1 - a2 R2 



where Nh = Vh = Uh = Ch = 0, while Bh is such that < Bh < l/a- 

Linearization of the equations about this fixed point leads to the following eigen- 
values of the matrix defining the linear part (adopting the order N, W, U, H, V, P, 
Q, k) 

70 = 1 (54) 



7i,2 = -^±^^ 



2 \ 4 1 - a^Bl 



11 



73,4 = -;±ji + aW (56) 




75 = (57) 

76 = -1 (58) 

77 = -2 (59) 

In particular, we see that the eigenvalues 71,2, which are related to the function W and 
determine the critical value a for the magnetic black holes [§, now depend on B^. The 
zero mode 75 is related to the occurrence of this free parameter. 

For Q = 0, the critical value of a corresponds to the transition of the eigenvalue 
71,2 (eq. ( |55D for Q = 0) from complex to real value Repeating this reasoning for 
dyonic black holes, we conclude, that the corresponding critical value of the coupling 
constant for the dyonic non-abelian black holes is given by 
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Here the relation of the critical value a to the charge Q seems only indirect, since 
the value of the function B{x) at the horizon enters in eq. (|60|) . However, for the dyonic 
non-abelian black holes a also represents the special point, where the non-abelian black 
hole solutions tend towards the corresponding extremal RN solution with horizon radius 
= ay/1 + Q"^ on the full interval x > x'f^ = xf^. This implies for the function B{x) 
according to eq. (|4^) 

lim,S(x) = — , yx>xr- (61) 

Consequently, we may insert for B^ in eq. (^) the constant RN value, eq. (|43|) , and 
obtain for the critical coupling constant the expression 



a{Q) 



3 + 402 

^ (62) 



which now depends on the electric charge Q. 

5.3 Fundamental dyonic black holes 

Let us now define some quantities helpful in the analysis of the numerical results. 
Denoting the set of functions of the extremal RN solution with horizon radius x^^ = 
ay/1 + Q'^ by (iV^N, -^^rn, -f^RN, -Prn), and denoting the set of functions of the non- 
abelian black hole solution with horizon radius Xh (with Xh < x]^^) by {N,K,H,P), 
we define 

DN = max^gr RN .|A/'RN(a;) - N{x)\ , (63) 
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and analogously we define the quantities DK, DP and DH. Clearly the quantities DN, 
DK, DP and DH must vanish, when the non-abelian black hole solutions bifurcate with 
the corresponding extremal RN solution. Therefore these quantities help to determine 
the critical values with good accuracy in the regime a > a. 
Q = 

Our numerical evaluation of the domain of existence of the magnetic non-abelian 
black holes in the a-Xh plane is in agreement with the results of and in particular 

we confirm a = v^/2. Considering the quantities DN, DK and DH (eq. (|63D) as 
functions of the horizon radius Xh, we observe that for a > a the three functions 
approach zero at a common value < a. As a — a changes sign, in particular the 

quantity DK changes drastically, becoming a convex function when the critical horizon 
radius x'f^ is approached. 

We have numerically determined the domain of existence of the fundamental dyonic 
non-abelian black hole solutions in the a-Xh plane, which is shown in Fig. 4 for Q = 1. 
The straight line indicates the extremal RN solutions, dividing the domain of existence 
into a lower region with non-abelian black holes only and an upper region with both 
non-abelian and non-extremal abelian black holes, completely analogously to the case 
of magnetic black holes. 

To start the detailed discussion of the diagram, we choose a small fixed value of 
the horizon radius Xh and vary a. Then the dyonic black hole solutions exist only up 
to some maximal value amax(a;/i). /,From amax{xh) a second branch extends backwards 
until a critical value adxh) is reached. With increasing Xh the second branch becomes 
increasingly smaller, and finally disappears. The critical values amax{xh) and adxh) 
are both shown in Fig. 4. 

The presence of the two branches is illustrated in Fig. 5 for Q = 1 and Xh = 0.8, 
where the quantities DN, DK, DP and DH are shown as functions of the parameter 
a. Here the second branch of solutions only exists in the small interval 1.1986 < a < 
1.2002. 

For fixed Xh and varying a, the functions of the black holes solutions evolve com- 
pletely analogously to those of the regular dyon solutions. Consequently, the limiting 
solution reached at the critical value adxh) consists of two parts, the outer part corre- 
sponding to the exterior of the extremal RN black hole with horizon xf;^ = ac\/l + Q^, 
unit magnetic charge and electric charge Q and the inner part representing a genuine 
non-abelian solution. 

Similarly for fixed a > a and varying Xh, the limiting solution reached consists 
of two parts, the outer part corresponding to the exterior of the extremal RN black 
hole with horizon x^^ = a\/Y^^f~Q^, unit magnetic charge and electric charge Q, as 
demonstrated in Figs. 6-7 for the black hole solutions with horizon Xh = 1, 1-1 and 
1.153, charge Q = 1 and a = 1 for the functions and P and the functions K and H, 
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respectively. 

According to eq. (|62D, for Q = 1 the critical value a is given by d = ^7/8 ~ 
0.935. Here the transition between the two regimes with different bifurcation patterns 
is supposed to occur. Indeed, our numerical calculations confirm this critical value, as 
demonstrated in Figs. 8-9, where the four quantities ^ DK , DP and DH are shown 
as functions of the horizon radius Xh for a = 0.94 and a = 0.93, respectively. 

For a = 0.94, which is slightly above a, the curves in Fig. 8 clearly approach 
zero at x'f^ ~ 1.235. This value is well below the horizon value of the corresponding 
extremal RN solution, x]^^ = a\/\ + ^ 1.329. Moreover, in this limit Xh — > x'^^^ the 
corresponding values of K{xh) and H{xh) do not tend to their RN values of zero and 
one, respectively, but they converge to generic numbers in the interval [0, 1]. 

In contrast, for a = 0.93, which is slightly below a, the function DK becomes a 
convex function when the critical horizon radius x'f^ ~ 0.93-\/2 ^ 1.315 is approached. 
Furthermore, the quantities K{xh) and H{xh) tend to their RN values zero and one, 
respectively, for Xh x'f^. Thus for a = 0.93 the limiting solution represents a RN 
solution on the full interval x > x". 

Clearly, the transition point a occurs inbetween these two values. Approaching a 
further from both sides, the value a = ^7/8 ~ 0.935 is confirmed. The gap in Fig. 4 on 
the rhs of a represents a tiny region, numerically not accessible with sufficient accuracy 
[P^. We have confirmed the validity of eq. (^) also for other values of the charge Q. 

In the regime a < a, we observe that the maximal value of the horizon radius Xh 
decreases with decreasing a, completely analogously to the magnetic case. Continuing 
our detailed discussion of the approach to the limiting solution in this regime, however, 
we observe differences to the magnetic case for smaller values of a and larger values of 
Q. 

For fixed a close to a and increasing Xh, the dyonic black hole solutions bifurcate 
with a RN solution when the maximal value of the horizon is reached. The limiting 
solution corresponds to the exterior of the RN solution on the full interval x > x". 
For smaller values of Q and somewhat smaller fixed a, with increasing Xh the solutions 
do not yet bifurcate with a non-extremal RN solution, when the maximal value of the 
horizon is reached. Instead a second branch of solutions emerges, extending backwards 
to lower values of the horizon radius, until it bifurcates with a non-extremal RN solution 
at the critical value of the horizon x'f^ |H, il Hi- Moving along both branches, the values 
of K{xh) and H{xh) change monotonically, reaching their respective RN values of zero 
and one as the critical horizon is reached. 

For the magnetic black holes as well as for the dyonic black holes with smaller values 
of Q this latter bifurcation pattern persists with decreasing a [§, If. In contrast, for 
the dyonic non-abelian black holes with larger values of Q new phenomena occur. 

We exhibit these new phenomena in Fig. 10, where the critical region of the a-Xh 
plane shown in Fig. 4 for Q = 1 is enlarged and more details are given. First, we 
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observe the occurrence of bifurcations, which lead to the small triangular region ABC, 
where three solutions exist. Thus below a we find four distinct regions in the a-Xh 
plane. With decreasing a and increasing Xh there are i) for 0.605 < a first one, then 
zero solutions, ii) for 0.575 < a < 0.605 first one, then three, then one, then zero 
solutions, iii) for 0.535 < a < 0.575 first one, then three, then two, then zero solutions 
and iv) for a < 0.535 first one, then two, then zero solutions. 

Second, with decreasing a, the value J^o of the limiting solution increases. When 
it reaches its maximal value of one, the function K{x) ceases to decay exponentially. 
The limiting solution then no longer represents a non-extremal RN solution, instead an 
oscillating solution is reached as the limiting solution. The corresponding critical point, 
where the transition from the limiting non-extremal RN solutions to the oscillating 
solutions occurs, is labelled D in Fig. 10. Along the curve BD the maximal value of 
the horizon radius is reached when the limiting value J^o = 1 is attained, whereas along 
the curves AB and A' A an intermediate extremal value of the horizon radius is reached 
when the limiting value J^o = 1 is attained. 

Figs. 11-13 illustrate both phenomena, representing the values K{xh), H{xh) and 
Joo as functions of the horizon radius for a = 0.7, 0.55 and 0.2, respectively. As seen 
in Fig. 11, for a = 0.7 only one branch of solutions exists. Here, with increasing Xh 
K{xh) decreases monotonically to zero, H{xh) increases monotonically to one, while Joo 
increases monotonically to a limiting value smaller that one. Decreasing a, the critical 
point C of Fig. 10 is reached, below which three branches of solutions exist. Moving 
along the three branches, K{xh) decreases monotonically, H{xh) increases monotoni- 
cally and Joo increases monotonically until it reaches one, as seen in Fig. 12 for a = 0.55. 
For still smaller values of a, at the critical value A, the third branch disappears, leav- 
ing two branches which end when J^o reaches one. This is seen in Fig. 13 for a = 0.2. 
The maximal horizon radius in this case is approximately 0.917. For small values of a 
(i.e. below the critical point A) the lower curve in Fig. 4 represents the critical value of 
the horizon, where the asymptotically exponentially decaying solutions cease to exist. 

Considering the Q-dependence of the new phenomena, we observe that the bifurca- 
tions leading to three branches arise only for Q > 0.8, and the limiting value Joo = 1 
is reached only for Q > 0.75. The latter is seen in Fig. 14, where the critical value D 
is shown, which marks the transition from the limiting non-extremal RN solutions to 
the oscillating solutions. Also shown in Fig. 14 is the critical value a. We note, that 
with increasing Q the two curves get increasingly closer. 

5.4 Excited dyonic black holes 

Besides the fundamental magnetic black hole solutions there exist radially excited black 
hole solutions, for which the gauge function K{x) possesses n nodes 0, ^. These 
radially excited solutions exist only for a < d = v^/2 [|, p. 
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Our numerical analysis of the radially excited dyonic black hole solutions strongly 
indicates that, at least for n = 1, the excited dyonic black hole solutions penetrate into 
a small domain of the a > a regime, as can already be anticipated from Fig. 3, where 
the first excited regular dyon solution is shown for charge Q = 1. Indeed, we observe 
that the maximal value of a, where the radially excited black hole solutions cease to 
exist, increases more strongly with Q than the critical value a. 

Let us now consider the domain of existence of the black hole solutions in the a — Xh 
plane in more detail. In the regime a > a, for fixed Xh and increasing a the solutions 
again reach a limiting solution at the critical value adxh), which consists of two parts, 
the outer part corresponding to the exterior of the extremal RN black hole with horizon 
x^^ = OcVl + Q^y unit magnetic charge and electric charge Q. 

For a — > ac{xh) the function N{x) develops a progressively decreasing minimum, 
A^min, reaching zero at x^^ = ac{xh)y/T+~Q^. This is illustrated in Fig. 15, where the 
value of A^min is shown as a function of a for the horizon radii Xh = 0.01, 0.4 and 0.8 
and charge Q = I. With increasing Xh the value of adxh) decreases, though it remains 
greater than a even for Xh = 0.8, as inspection of the four quantities DN, DK, DP 
and DH reveals. 

In Fig. 16 we show the functions N{x), K{x), H{x) and P{x) of the first radially 
excited black hole solution for charge Q = 1, horizon radius Xh = 0.4 and a = 0.938 > a. 
The minimum of the function N{x), which here occurs close to x]^^ = ac{xh)V^ + Q^: 
remains well separated from the horizon Xh for a adxh)- 

Considering finally the first excited black hole solutions for small fixed a and in- 
creasing Xh, we observe that as for the fundamental black hole solutions, the quantity 
Joo here plays a major role. With increasing Xh, Joo approaches its maximal value one, 
where the asymptotically exponentially decaying solutions cease to exist. For a = 0.2, 
for instance, this happens at Xh ~ 0.28. 

6 Conclusions 

In analogy to gravitating monopoles ^ also gravitating dyons exist. When 
a regular horizon Xh is imposed, from these solutions "black holes within monopoles" 
[^, ^ and "black holes within dyons" emerge. Besides the fundamental regular 
and black hole solutions also radially excited solutions exist [@, |], related to the 
corresponding EYM solutions [|14|, |1^. The static spherically symmetric "black holes 
within monopoles" and "black holes within dyons" provide counterexamples to the 
"no-hair conjecture". 

The domain of existence of the fundamental dyonic non-abelian black hole solutions 
in the a — Xh plane is similar to the one of the fundamental magnetic non-abelian black 
hole solutions. There are two regimes of the coupling constant a. For a < a, the 
maximal horizon radius increases with increasing a, whereas for a > a it decreases. 
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The transition point between these two regimes, d, depends on the charge Q of the 



solutions, d = ^(3 + AQ^)/{A{1 + Q2)). 

For fixed a and varying horizon radius, a hmiting solution is approached. For a > a 
the limiting solution consists of two distinct parts. The outer part corresponds to the 
exterior of an extremal RN black hole solution with horizon radius x]^^ = a\/\ + Q^^ 
unit magnetic charge and electric charge Q, whereas on the interval xy,'^x < x^ the 
limiting solution has the features of generic non-abelian solutions. At the transition 
point a = a the limiting solution corresponds precisely to the exterior of the extremal 
RN black hole solution with horizon radius x^ = ay/1 + Q'^, on the full interval 
X > Xh- In contrast, for a < a the limiting solution represents the exterior of a non- 
extremal RN solution for the larger values of a, while for the smaller values of a the 
dyonic non-abelian black hole solutions encounter a critical point, where the solutions 
no longer are exponentially decaying but become oscillating. 

The domain of existence of the radially excited dyonic non-abelian black hole so- 
lutions in the a — Xh plane is similar to the one of the fundamental solutions. Again, 
there are two regimes of the coupling constant a with the transition point a. This is in 
contrast to the radially excited magnetic non-abelian black hole solutions, which exist 
only up to q; 0, |[. 

In a large part of their domain of existence, the magnetic non-abelian black hole 
solutions are classically stable 0, |^, H, |19[. In contrast, dyonic non-abelian black hole 
solutions should be classically unstable. As for the classical dyon solutions in fiat space, 
their mass should be lowered continuously, by lowering the electric charge, as long as 
there is no charge quantization p|. 
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7 Appendix A 



We show the local existence of the non-abelian dyon solutions at the singular point 
X = in Appendix A.l and the local existence of the non-abelian dyonic black hole 
solutions at the singular point Xh in Appendix A. 2. 



7.1 Appendix A.l 

As for the monopoles, the field equations for the dyons have singular points at a; = 
and X = oo. For the monopoles there exists a two-parameter family of local solutions 
regular at x = 0. Analogously, for the dyons there exists a three-parameter family of 
local solutions regular at x = 0, as shown in the following. 

Let us consider eqs. ([T7|)-(|?I|), using notation I. At the singular point x = 0, the 
expansion of regular solutions depends on four parameters, yielding 

N{x) = 1 - nax^ + O(x^) , 

A{x) = Ao{l + a2x'^ + 0{x^)) , 

K{x) = 1 - A;2X^ + O(x^) , (64) 

J(x) = Aoijix + 0{x^)) , 

H{x) = hix + O(x^) , 

where Aq, hi, ji and k2 are arbitrary parameters while 

n2 = a\jl + hl + Akl + ^) . (65) 

6 

Eliminating the metric function A, we now show that the remaining equations for 
A^, K, H and P admit a three-parameter family of local solutions near x = 0, analytic 
in the parameters. This is the counterpart of Proposition 2 of Appendix A of @ in the 
presence of the dyon degree of freedom. Let us define 



H l-K NK' 

Wi = — , W2 = NH , W3 = — , W4 = , 

P I - N 

w, = - , We = iP' + PM)N , = — (66) 

X x^ 

with 
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and furthermore 

ui = ^{2wi + W2) , U2 = ^{ws-Wi), U3 = ^{2w5 + We) , 

^1 = ^(^1-^2), f2 = ^(2^3 + ^4) , V3 = ^{w5-Wq), 

^4 = ^ - ^(«? - ^vl + Aul - 2vl + ul- 2vl + ^) . (68) 
Then eqs. (^) and (|TPp-(^Tp can be rewritten in the form 

CCUj^ — X ^ 1 — X 3 ^ 

xv'i = -XiVi + x^Qi , i = 1 - 4 , (69) 

where fi and gj are analytical functions of and and A, = 3, z = 1, 2, 3. It 

now follows from Proposition 1 of Appendix A of that the system of equations ([T7|) 
and (|l^)-(^) admits a three-parameter family of solutions of the form 

Ui = Ci + 0{x^) , V, = 0{x^) (70) 

which are locally analytic in and in the constants q. These constants correspond to 
hi, ji and ^2 in (|6^) . 



7.2 Appendix A. 2 

Following [H and considering the system of equations (^7D -(|35D, we here demonstrate the 
local existence of the non-abelian dyonic black hole solutions at the singular point Xh- 
Black hole solutions are characterized by a first order pole singularity of the function k 
at the non-degenerate horizon p[. Due to translation invariance of the equations under 
T —>■ T + tq, we can assume that the horizon occurs at r = 0. 

To rewrite the set of equations in the form (|69D, which reveals the number of free 
parameters occurring in the expansion at the horizon of the locally analytic solution, 
we introduce the following set of functions 

(Ml, U2, Us, M4, Ms) = (r, W, H, 2kN - 2U^ - V\ S) , 

{vi, V2, V3, v^) = {k--,U,V,Y- - -)) , (71) 

r r 

with 

S = B/t , Y = S . (72) 

These functions are regular at the horizon and can be expanded in powers of r. 
Moreover, it follows that the equations for Ui and Vk are of the form (|69D with Ai = 2, 
A2 = A3 = 1 and A4 = 3. This demonstrates the existence of a five-parameter family 
of local analytic solutions of the equations. Because of the constraint (|36|), it actually 
reduces to a four-parameter family of local solutions. 
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Figure Captions 
Figure 1 

The normalized mass of the regular dyon solutions is shown as a function of the 
charge Q for vanishing Higgs self-coupling {13 = 0) and finite Higgs self-coupling = 
0.5) for the fiat space solutions (solid lines) and curved space solutions with a = 0.5 
(dashed lines). The asterisks mark the transition from exponentially decaying solutions 
to oscillating solutions. 
Figure 2 

Same as Fig. 1 for the asymptotic value of the function J{x). 
Figure 3 

The mass of the solutions of the 1st excited dyon branch for Q = 1 and the 1st 
excited monopole branch {Q = 0) is shown as a function of the coupling constant a for 
P — (sohd) . Also shown is the mass of the branch of extremal RN solutions with unit 
magnetic charge and Q = 1 as well as Q = (dotted). The normalization is chosen 
such that at q; = the mass of the first Bartnik-McKinnon solution is obtained. 
Figure 4 

The domain of existence of the dyonic black hole solutions in the a — Xh plane is 
shown for Q = 1. The straight diagonal line represents the extremal RN solutions. To 
the right of this line only non-abelian black hole solutions exist. Here the limiting curve 
represents the maximal values of a, while the second curve close to it represents the 
critical values of a, where the solutions bifurcate with the corresponding extremal RN 
solution. The asterisk marks the critical value a. To the left of the line of extremal RN 
solutions both non-extremal RN solutions and non-abelian black hole solutions coexist 
below the upper curve, representing the maximal horizon radius as a function of a. 
The small rectangle indicates the critical region, where bifurcations occur (enlarged in 
Fig. 10). Below the bifurcations, the lower curve represents the critical values of the 
horizon, where the solutions terminate because Joo = 1 is reached. 
Figure 5 

The quantities DN, DK, DP and DH are shown as functions of the parameter a 
for the dyonic black hole solutions with Q = 1 and Xh = 0.8. The second branch of 
solutions in the interval 1.1986 < a < 1.2002 is clearly visible. 
Figure 6 

The functions N{x) and P{x) are shown for the dyonic black hole solutions with 
Q — 1 and q; = 1 for the horizon radii x^ — 1, 1.1 and 1.153. 
Figure 7 

Same as Fig. 6 for the functions K[x) and H[x). 
Figure 8 

The quantities DN, DK, DP and DH are shown as functions of the horizon radius 
Xh for the dyonic black hole solutions with Q —1 and a — 0.94 > a. 
Figure 9 
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Same as Fig. 8 for a = 0.93 < a. 
Figure 10 

Enlargement of the critical region (small rectangle) of Fig. 4, showing the domain 
of existence of the dyonic black hole solutions in the a — Xh plane for Q = 1. The 
numbers 0-3 indicate the numbers of solutions in the respective areas. A and C indicate 
the critical values, where the bifurcations occur, D indicates the critical value, where 
the transition from the limiting non-extremal RN solutions to the oscillating solutions 
occurs. Joo = 1 is reached along the curve A'ABD. 
Figure 11 

The quantities K{xh), H{xh) and Jqo are shown as functions of the horizon radius 
Xh for the dyonic black hole solutions with Q = 1 and o; = 0.7. 
Figure 12 

Same as Fig. 11 for a — 0.55. 
Figure 13 

Same as Fig. 11 for a — 0.2. 
Figure 14 

The charge Q, corresponding to the critical value D, where the transition from the 
limiting non-extremal RN solutions to the oscillating solutions occurs, is shown as a 
function of a (dotted). Also shown is the critical value a (sohd). 
Figure 15 

The minimum of the function N[x) is shown as a function of the parameter a for 
the first radially excited dyonic black hole solutions with Q = 1 and horizon radii 
Xh = 0.01, 0.4 and 0.8. 
Figure 16 

The functions N{x), P{x), K{x) and H{x) arc shown for the first radially excited 
dyonic black hole solution with Q = 1, a = 0.938 and Xh = 0.4. 
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